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Abstract 

The cross sections for e~^e~ (j)iT~^Tr~ and e~^e~ 0/o(98O) arc measured from threshold to 
^/s = 3.0 GeV using initial state radiation. The analysis is based on a data sample of 673 fb~^ 
collected on and below the T(45') resonance with the Belle detector at the KEKB asymmetric- 
energy e'^e~ collider. First measurements are reported for the resonance parameters of the ^(1680) 
in the ^tt+tt" mode: m = (1689 ± 7 ± 10) McV/c^ and T = (211 ± 14 ± 19) MeV/cP. A structure 
at ^/s = 2.1 GcV/c^, corresponding to the so called y(2175), is observed; its mass and width are 
determined to be 2079 ± 13l^| MeV/c^ and 192 ± 23l^f MeV/c^, respectively. 

PACS numbers: 13.66.Bc, 13.25.-k, 14.40.Cs 
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Although vector mesons are produced copiously in e"'"e~ annihilation, the resonance 
parameters of excited vector states are not well measured [l|. In the ss sector, the 
0(1680) was first observed by the DM1 Collaboration in the reaction e"'"e~ KKir ^ 
and was recently studied in the initial-state radiation (ISR) events of the e~^e~ (prj and 
K* {^92)K + c.c. modes by the BaBar Collaboration, which measured its mass and width to 
be 1723 ± 20 MeV/c^ and 371 ± 75 MeV/c^, respectively (s"]. In a study of ISR events of 
the type, e^e~ TisR^^^^i""? the BaBar Collaboration observed two clear structures near 
^/s = 1.68 GeV/c^ and 2.175 GeV/c^, the latter was produced dominantly via a 0/o(98O) 
intermediate state, and was dubbed the F(2175) fl. The BES Collaboration confirmed the 
y(2175) in the 0/o(98O) invariant mass spectrum of J/ip ?70/o(98O) [r] —>■ 77, (p —>■ K^K~ 
and /o(980) — > tt+tt") decays with a statistical significance of about 5(T [5|. The Particle 
Data Group (PDG) assigns all these observations to a new state referred to as the 0(2170) [l[. 

Since the y(2175) resonance is produced via ISR in e+e~ collisions, its J^*" = 1 . 
This observation stimulated the theoretical speculation that y(2175) may be an s-quark 
counterpart of the F(4260) 0, since both are produced in e^e annihilation and exhibit 
similar decay patterns. On the other hand, a number of different interpretations have been 
proposed for the y(2175) with predicted masses that are consistent, within errorSjWith the 
experimental measurements. These include: an ssg hybrid jl]; a 2^Di ss state with a 
width predicted to be in the range 120-210 MeV/c^; a tetraquark state llol,|ll|; a AA bound 
state [l2|; or an ordinary 0/o(98O) resonance produced by interactions between the final 
state particles [l3|. The possibility that the y(2175) is a ?)^Si ss state is disfavored by the 
rather large predicted width (F ~ 380 MeV/c^) |l4l|. which disagrees with the experimental 
upper limit of 100 MeV/c^ [4]. A recent review [l5| discusses the basic problem of the large 
expected decay widths into two mesons, which contrasts with experimental observations. 

In the analysis reported here, we use a data sample with an integrated luminosity of 
673 fb~^ collected with the Belle detector [3] operating at the KEKB asymmetric-energy 
e+e^ (3.5 on 8 GeV) collider [13] to investigate the (pix^Ti" final state via ISR. About 90% 
of data were collected at the T(4S') resonance (a/s = 10.58 GeV), and the rest were taken 
at a center-of-mass (CM.) energy that is 60 MeV below the T(4S') peak. 

For Monte Carlo (MC) simulation of the ISR process, we generate signal events with 
the PHOKHARA program In this program, after one or two photons are emitted, the 
lower energy e^e~ pair forms a resonance X that subsequently decays to ^tt+tt" or 0/o(98O) 
with the (/> decaying into K^K~ and the /o(980) decaying into tt+tt^. In the X —>■ ^tt+tt" 
generation, we assume that the tttc system is pure S-wave and that it and the (p are also in 
a relative S-wave. The vr+vr" invariant mass distributions are generated according to phase 
space. 

To select the ^vr^vr^ final states, we use the following event selection criteria. We require 
four good charged tracks with net charge equal to zero. Good tracks should have impact 
parameters perpendicular to and along the beam direction with respect to the interaction 
point of less than 0.3 and 2 cm, respectively, and transverse momentum is restricted to be 
higher than 0.1 GeV/c. For each charged track, the particle identification information from 
different detector subsystems is combined to form a likelihood for each particle species (i), 
Li 1^. Tracks with 71^ = lI^'^l < 0.9 are identified as pions, tracks with TZk > 0.6 



are identified as kaons. We require that two tracks be identified as pions, and assign kaon 
masses to the other two tracks, of which only one is required to be identified as a kaon. The 
sum of the charged track energies and that of reconstructed photons (neutral clusters in the 
electromagnetic calorimeter with energies greater than 100 MeV and not associated with 
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charged tracks in the central drift chamber) must be greater than 9 GeV; this ensures that an 
ISR photon is detected. The efficiency of this requirement is nearly 100% for 07? "•"vr" invariant 
masses below 3.0 GeV/c^. We require -1.0 (GeV/c^)^ < M"^^^ < 2.0 (GeV/c^)^ to remove 
multi-photon backgrounds, where M^g^ is the square of the mass that is recoiling against the 
four charged tracks. The MC simulation indicates that the application of these requirements 
reduces the contamination from the background ISR process e^e~' — ^ 07r"'"7r~7r°(7) jlO] to 
the 1% level; the effect of the residual background is included as a systematic error. 

In events with K'tt^h' mass below 3.0 GeV/c^, a (f) meson signal is evident in the 
K^K~ invariant mass distribution, as shown in Fig. [H This distribution is fitted with a 
Breit-Wigner (BW) (the meson mass and width fixed to PDG values) convolved with 
a Gaussian resolution function as the signal and a second-order polynomial as the back- 
ground shape; here the width of the Gaussian resolution function is fixed at a = 0.95 MeV/c^, 
its MC-determined value. The fit yields 4832 ± 132 (pir'^ir' events. 

1200^ 
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FIG. 1: The K^K~ invariant mass distribution for selected e'^e~ K^ir^ir^ candidates 

with K^ir^TT^ invariant mass less than 3.0 GeV/c^. The curves show the fit projections for 
the full fit and for the background component. 

We obtain the number of 07r+7r^ events in each K'^K~tt^tt~' invariant mass bin by fitting 
the K^K~ invariant mass distribution for each interval, using the same fit model as in the 
fit for the overall sample. In order to reliably control the background shape, we restrict the 
coefficients of the background polynomials in nearby bins to vary smoothly along parabolas. 
The parameters of these parabolas are determined from fits to the coefficients obtained from 
fits to the K'^K~ invariant mass distribution in each K^K~7t'^it~ invariant mass bin. The 
resulting 07r+7r~ invariant mass distribution is shown in Fig.[2]^a). Combinatorial background 
has a smooth distribution in the K^K~ invariant mass distribution and does not affect the 
results since we fit the K^K~ distribution to obtain the number of (j)7T^7r~ events in each 
K^K~7^~^7r~ mass bin. In Fig. [2]^a) there are two distinct peaks: one near 1.7 GeV/c^ and 
another near 2.1 GeV/c^, corresponding to the 0(1680) and y(2175) states. In addition, 
there is a cluster of events near 2.4 GeV/c^. 
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FIG. 2: The ^vr^vr^ invariant mass distribution obtained from fitting to ra{K^ K^) in each 
m{K~^K^TT^Tr^) bin (a) as described in the text, and the measured e~^e~ — > (/ivr^vr" cross sec- 
tion (b). The errors are statistical only. 



The e^e (p-n^-n cross section for each ^tt^tt mass bin is computed using 



where nf^, Si, and Ci are the number of 07r+-7r~ events fitted in data, the efficiency, and 
the effective luminosity j2l| in the i-th ^vr+vr" mass bin, respectively; B{(j) K^K~) = 
49.2% According to the MC simulation, the efficiency increases smoothly from 1.72% to 
2.67% for (f)TT~^TT~ masses ranging from threshold to 3.0 GeV/c^. The resulting cross sections 
are shown in Fig. [2](b) [22*] . 

We now consider the quasi-two-body intermediate state 0/o(98O). Figure [31(a) shows a 
scatter plot of m{(f)7i~^7r~) versus m{n^7r~) for events in the (p signal region (1.013 GeV/c^ < 
iTT'K+K- < 1-025 GeV/c^). Non-0 background is subtracted according to the normalized (p 
mass sidebands {niK+K- e [0.989, 1.013] GeV/c^ or ttik+k- e [1-025, 1.049] GeV/c^). It can 
be seen that within the /o(980) mass band there are clusters of events that correspond to 
F(2175) ^ (^/o(980) and J/ip 0/o(98O). Figure[3](b) shows the projected m(7r+7r-) distri- 
bution, where a clear /o(980) signal is evident. For each 25 MeV/c^ bin of m{K'^ K~7r~^7r~) , 
we select events with m(7r+7r^) G [0.85, 1.10] GeV/c^ [23], and fit their m{K^ K^) distribu- 
tion to extract the number of 0/o(98O) events in a way similar to that used to extract the 
number of (j)TT^TT~ events described above. The resulting (/)/o(980) invariant mass distribu- 
tion is shown in Fig. Il](a). A very clear y(2175) signal is observed with an accumulation of 
events around 2.4 GeV/c^. The e'^e~ — > 0/o(98O) cross section for each 0/o(98O) mass bin 
is computed in the same way as that for the 07r+7r^ mode using i3(/o(980) tt+tt^) = |. 
The efficiency is nearly independent of 0/o(98O) mass, 2.3% from threshold to 3.0 GeV/c . 
The resulting cross sections are shown in Fig. 111(b). 

The sources of the systematic errors for the cross section measurements are summarized in 
Table [H The uncertainty is 1.0% for pion identification, and is negligible for kaons since the 
identification of only one of the kaons is required; the uncertainty in the tracking efficiency 
is 1% per track, and is additive; the efficiency uncertainty associated with the M"^^^ require- 
ment is determined to be 1.0% [3]; the uncertainties in the yields of (/)7r"'"7r~ and 0/o(98O) 
events for each mass bin due to the signal fit are estimated to be 5.0% and 4.0% as deter- 
mined by: changing the mass resolution, the orders of the background polynomial, and the 
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FIG. 3: Scatter plot of m{4>7r~^TT~) versus m(7r+7r^), where non-cp background is subtracted ac- 
cording to the normahzed 4> mass sidebands (a), and the projection on 7r"'"7r~ mass, a clear /o(980) 
signal is visible (b). 
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FIG. 4: The 0/o(98O) invariant mass distribution (a) obtained from fitting to m{K^ K^) in each 
m{K~^ K~ fQ{9SQ)) bin as described in the text, and (b) the measured e+e" — > (/>/o(980) cross 
sections. The errors are statistical only. 



fitting range. The error due to the uncertainty of the background parameters changes with 
m^TT+TT- (GeV/c^) because of signal-to-noise ratio variations between threshold and 3 GeV/c^. 
The error is approximated as (5.92 — 5.44m(^7r+7r- + l-28m^^+^_)/(0.02 + (m^^+Tr- — 1.9)^)%, 
m<^^+^- in GeV/c^, which is about 4% in the 0(1680) and the F(2175) mass regions, and in- 
creases continuously to about 10% in between and decreases to about 1% above 2.3 GeV/c^. 
This is common to the 07r+7r~ and (/)/o(980) modes. The </)7r+7r^7r°(7) background fur- 
ther contributes a 1% uncertainty to the (j)7i^Tr~ , and a negligible fraction to the 0/o(98O) 
mode. Belle measures luminosity with a precision of 1.4% using wide angle Bhabha events. 



The uncertainty of the ISR photon radiator is 0.1% [2l|. The main uncertainty in the 



PHOKHARA [l8j generator is due to the modelling of the tt+tt^ mass spectrum. This is 
tested by generating events with different tt+tt^ mass distributions. We take 5% and 3% 
as conservative uncertainties related to the MC generators for the 07r+-7r~ and 0/o(98O) 
modes, respectively. According to the MC simulation, the trigger efficiency for these events 
is around (97 ± 1)% with little dependence on the ^tt+tt" invariant mass. The uncertainty 
of B{(p K+K-) is 1.2% Finally the error due to MC statistics is 0.8%. We assume all 
the sources are independent and add them (excluding that from background parameters) in 
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quadrature, resulting in the total systematic errors on the cross sections of 8.6% and 6.9% 
for 07r+7r~ and 0/o(98O) modes, respectively. 

TABLE I: Systematic errors (%) in the (pn^n^ and (j)fo(980) cross sections measurements. 
The systematic error from the uncertainty of background parameters [(5.92 — 5.44m^^+^- + 
1.28m2 +^_)/(0.02 + (m^^+^- - 1.9)2)%, m^^+^- in GeV/c^] is not included. 



Source 


(pTT^TT 


</./o(980) 


Particle ID 


1.0 


1.0 


Tracking 


4.0 


4.0 


M^^gj, selection 


1.0 


1.0 


(p signal fit 


5.0 


4.0 


(pTT^ TT^ TT^ ('j) background 


1.0 


0.0 


Integrated luminosity 


1.4 


1.4 


Generator 


5.0 


3.0 


Trigger efficiency 


1.0 


1.0 


Branching fractions 


1.2 


1.2 


MC statistics 


0.8 


0.8 


Sum in quadrature 


8.6 


6.9 



There are strong J/ip signals in both 07r+7r~ and 0/o(98O) samples (see Fig. [5]). 254 ±23 
J/ip ^ (f)7i^7i^ and 60 ± 11 J/ip (/)/o(980) events are obtained from fitting the J/ijj 
signals and subtracting the peaking backgrounds by fitting the normalized sidebands 
events. We determine B{J/ip ^Tr+Tr") x T^+e- = (4.50 ±0.41 ±0.26) eV/c^ and i3(J/V' ^ 
0/o(98O)) X Te+e- = (2.22 ±0.41 ±0.13) eV/c^, where the first errors are statistical, and the 
second systematic. The sources of the systematic errors are almost the same as those listed 
in Table [H except that the uncertainty of the signal fit is replaced by that of the JM 
signal fit. Our results are consistent with the PDG values pj] and BaBar's measurements 
but are more precise. 




2.9 3 3.1 3.2 3.3 2.9 3 3.1 3.2 3.3 



m((t) n n) (GeV/c^) m(0 fg) (GeV/c^) 

FIG. 5: The J/ip — > ^vr+vr^ (a) and J/ip (/)/o(980) (b) signals in the (p signals regions (open 
histograms) and in the 4> sidebands (shaded histograms, normalized). The curves are the best fits 
to the signal and sidebands distributions, respectively. 
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In order to obtain the resonance parameters of the 0(1680) and y(2175), a least squares 
fit is apphed to the ^tt+tt" cross section distribution. Since the 0(1680) decays into 071"*" vr" 
while the F(2175) decays dominantly into 0/o(98O), we use two incoherent BW functions 
in the fit, one for the 0(1680) and the other for the y(2175). The fit result is shown 
in Fig. El^a), with a goodness-of-the-fit of I'^^df = 68/55, corresponding to a C.L. of 
12%. The statistical significance of each resonance is greater than lOcr. From the fit we 
obtain the following resonance parameters of the 0(1680): m = (1689 ± 7) MeV/c^, F = 
(211 ± 14) MeV/c^ and B{^7i+7r-) x Be+e- = (1-24 ± 0.09) x 10"^ while those of the 
r(2175) are m = (2079 ± 13) MeV/c^, F = (192 ± 23) MeV/c^ and B{(pn+n-) x Be+e- = 
(1.10 ± 0.10) X 10~^, where the errors are statistical only. We also perform a fit with an 
additional incoherent BW function centered near 2.4 GeV/c^. The fitted parameters of this 
structure are m = (2406 ±32) MeV/c^ and F = (57 ±58) MeV/c^ with a goodness-of-the-fit 
of /ndf = 62/52, corresponding to a C.L. of 15%, where the errors are statistical only. 
The statistical significance of the structure at 2.4 GeV/c^ is 1.5cr as determined from the 
change in the value. A fit with an additional non-resonant component does not improve 
the fit quality, and the contribution of the non-resonant term is negligibly small. 

We fit the 0/o(98O) cross section distribution with a single BW function that interferes 
with a non- resonant component which is partly from non-0/o(98O) contribution, and partly 
from the possible 0(1680) — > 0/o(98O) at the high mass tail of the 0(1680), as in BaBar's 
analysis There are two solutions with very similar resonance parameters. The interfer- 
ence is constructive for one solution and destructive for the other. Figure [6](b) shows the 
result. The fit yields m = (2163±32) MeV/c^ and F = (125±40) MeV/c^, with a goodness- 
of-the-fit of x'^/^^df = 49/36, corresponding to a C.L. of 8.0%. Here we quote simple averages 
of the two solutions, and enlarge the errors to cover the full uncertainties. We also perform 
a fit with an additional coherent BW function near 2.4 GeV/c^. The fitted mass and width 
of the y(2175) are consistent with the results listed above, and the statistical significance 
of the structure at 2.4 GeV/c^ is estimated to be 2.3a. 
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FIG. 6: Fit to (a) e^e~ — > 07r^7r^ cross section with two incoherent BW functions, one for the 
0(1680) and the other for the y(2175) and (b) e^e~ 0/o(98O) cross section with a single BW 
function that interferes with a non-resonant component. The curves show the projections from the 
best fit and the contribution from each component. In (b), the dashed curves are for the destructive 
interference solution and the dot-dashed curves for the constructive interference solution 



Since the differences in the y(2175) resonance parameters from fits to 07r+7r and 
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(/>/o(980) are due to the assumptions on the background shape and the existence of ad- 
ditional resonances, and 0/o(98O) events are a subsample of ^vr+vr" events, we take the 
values of the y(2175) resonance parameters from the fit to the ^vr+vr" cross section with 
two incoherent resonances as the central values and the differences are taken as one source 
of systematic errors, which is the dominant one for the y(2175) resonance parameters. For 
the 0(1680) and y(2175) resonance parameters, we have also considered the uncertainties 
in the absolute mass scale, the mass resolution, the parameterization of the resonances 
and background shape, fit range and possible existence of additional resonances as sys- 
tematic errors. Finally, we obtain m(0(168O)) = (1689 ± 7 ± 10) MeV/c^, r(0(168O)) = 
(211 ± 14 ± 19) MeV/c2, and S(0(168O) ^ ^tt+tt") x B^+e- = (1-24 ± 0.09 ± 0.14) x 10"^ 
and m(r(2175)) = (2079 ± IStH) MeV/c^, r(F(2175)) = (192 ± 23^11) MeV/c^, and 
;B(F(2175) 07r+7r-) x B^+e- = (1-10 ± 0.10 ± 0.12) x 10"^ 0], where the first errors are 
statistical, the second systematic. 

In summary, we present the most precise measurements of the cross sections for e~^e~ — >■ 
07r+7r^ and e~^e~ 0/o(98O) from threshold to y/s = 3.0 GeV. The masses and widths 
of the 0(1680) and y(2175) are determined and are in agreement with the previous mea- 
surements [H, 0, Si, S- The width of the F(2175) tends to be lar ger than in previous 
measurements [4, l5| although the error is large. We find that the widths of the 0(1680) 
and y(2175) are quite similar and both are at the 200 MeV/c^ level. This may suggest 
that the y(2175) is an excited 1 ss state. Since the /o(980) is thought to have a large 
ss component, F(2175) 0/o(98O) can be viewed as an open-flavor decay as opposed to 
the case of y(4260) J/tp7T~^Tr~ , which is a hadronic transition. The study of the y(2175) 
in other decay modes would be useful for distinguishing between different possibilities. The 
branching fractions of J/ip decays into 07r+7r~ and 0/o(98O) are measured with improved 
precision; the results are in good agreement with the existing results. 
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